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The structural and electronic properties of cubic GaN are studied within the local 
density approximation by the full-potential linear muffin-tin orbitals method. The 



the potential and charge density. The influence of d electrons on the band structure, 



Ga 3d electrons are treated as band states, and no shape approximation is made to 
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■ charge density, and bonding properties is analyzed. It is found that due to the 

<n : 

■ energy resonance of the Ga 3d states with nitrogen 2s states, the cation d bands 



are not inert, and features unusual for a III-V compound are found in the lower 



■ part of the valence band and in the valence charge density and density of states. 

To clarify the influence of the Ga d states on the cohesive properties, additional full 
and frozen-overlapped-core calculations were performed for GaN, cubic ZnS, GaAs, 

x 

5_j ■ and Si. The results show, in addition to the known importance of non-linear core- 

valence exchange-correlation corrections, that an explicit description of closed-shell 
repulsion effects is necessary to obtain accurate results for GaN and similar systems. 
In summary, GaN appears to be somewhat exceptional among the III-V compounds 
and reminiscent of II- VI materials, in that its band structure and cohesive properties 
are sensitive to a proper treatment of the cation d bands, as a result of the presence 
of the latter in the valence band range. 
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I. INTRODUCTION 



Recently, considerable interest has arisen |IJ in the wide-gap III-V nitrides as candidates for 
the realization of light emitting devices in the blue range of the visible spectrum. In particu- 
lar, succesful epitaxial growth of thin films of gallium nitride has recently been demonstrated, 
resulting in materials having the wurtzite or the zincblende structure ||, depending on the 
substrate on which the growth was effected. 

Despite its increasing technological interest, GaN is poorly known from the theoretical point of 
view. Accurate theoretical predictions are of some relevance in this case, since characterization 
and experimental work on this material is still at an early stage when compared to the average 
developement level of III-V technology. The issue of a theoretical description of this system 
is also of interest in itself, due to the open question about the role of the Ga d electrons. 
In this work we present calculations of lattice properties and band structure of GaN in the 
zincblende structure (ZB) performed ab initio within the local-density approximation (LDA) 
to density functional theory [[| using the all-electron full-potential linear muffin-tin orbitals 
method (FP-LMTO) [[|. We analyze the character of band states, charge densities and density 
of states, showing the relevance of the Ga d shell in determining the properties of the material. 
To demonstrate the effects of the d shell on cohesive properties, we perform full and frozen- 
overlapped-core calculations for GaN, ZnS, GaAs, and Si monitoring the role of the d electrons 
and of core relaxation in the four materials. 

In Sec. II we discuss technicalities, in Sec. Ill some background, in Sec. IV the results are pre- 
sented and discussed, and in Sec.V we examine the effects of d electrons on cohesive properties, 
charge density and density of states (DOS). Rydberg and bohr atomic units are used throughout 
the paper unless otherwise stated. 

II. TECHNICALITIES 

Our calculations were performed within the LDA by the FP-LMTO method |4j using the 
Vosko-Wilk-Nusair || parametrization of the Ceperley-Alder || exchange correlation energy, 
and the Monkhorst-Pack 10 special-points mesh J7j. The 2s and 2p states of N and the 3d, 4s, 



and 4p states of Ga are treated as bands, whereas the remaining core states are self-consistently 
relaxed in a spherical approximation. As variational basis for solving the Schrodinger equation, 
three augmented Hankel functions are used on each occupied atomic site with decay energies 
-0.7, -1.0 and -2.3 Ryd up to angular momentum I = 2. This LMTO basis is centered only on 
the Ga and N spheres, giving 27 functions per atom. 

The atomic spheres are non-overlapping (sphere radii are 97 % of half the nearest-neighbor 
distance) in contrast with the often-used atomic-spheres approximation (ASA) version of the 
LMTO method §|. As usual, empty spheres (of the same size as the atomic ones) are inserted 
in the interstitial regions of the ZB structure to improve the packing fraction. 

A full-potential technique needs an accurate evaluation of interstitial three-center integrals 
and charge density. This is done by an interpolation technique which represents the product 
of two Hankel functions as a linear combination of other Hankel functions in the interstitial 
region. The same technique is used to represent the interstitial exchange-correlation potential 
and energy density. More details are to be found in Ref. [|J. The auxiliary interstitial charge 
density is expressed as an expansion in Hankel functions with decay energies -1.0 and -3.0 Ryd, 
centered in all spheres, with angular momentum up to / = 4. 

The result presented below were obtained using a non-relativistic code. We also performed a 
scalar-relativistic calculation for GaN with the same method, and the results are found to be 
essentially unchanged (the lattice constant and fundamental gap are reduced by 0.2 % and 1.5 
% respectively, and the cohesive energy is practically identical). 



III. BACKGROUND 



Although GaN is close to being a polytypic material |],|9| , we have not attempted a prediction 



of the relative stability of ZB and wurtzite structures |10[. Such a prediction would be somewhat 
secondary to our understanding of this system. The first reason for this is that wurtzite-ZB 
energy differences are known to be very small both experimentally and computationally in 



many analogous systems |TTJ; second, wurtzite is already known to be the stable equilibrium 
structure of GaN under normal experimental conditions; finally, the existence of an "easy" 
structural transitions induced by the substrate symmetry in epitaxial growth has already been 



demonstrated experimentally |],|2|]. Therefore, and also in view of the existence of a synthetic 
cubic phase of GaN, we hold an analysis of the properties of the cubic phase for more relevant 
at this stage. 

In view of the real-space localization of the valence electrons of N and of the d shell of Ga, the 
use of an all-electron method, in which electronic states are all treated on equal footing irrespec- 
tive of their localization properties, appears to be important. In particular, prior to information 
to the contrary, a proper treatment of the Ga d shell is required. To get an indication of the 
relevance of these states, one can consider calculations for the free atoms. The approximation 
of freezing any electronic states in the core, or of eliminating them by a pseudization procedure, 
is only plausible if they are clearly more strongly bound than any of the relevant valence states 
in the free atom. Further, even if the frozen-core or pseudoatom approximation appears to 
be reasonable from inspection of the free atoms, it may break down in a compound crystal in 
the presence of resonances of the atomic levels of different constituent atoms. An inspection 
of the free-atom energy levels of Ga and N shows that a special situation occurs in GaN (see 
Fig.|I|). The upper valence bands will originate from the sp states of Ga and the 2p electrons 
of N, the corresponding non-relativistic full-core free-atom eigenvalues being e£f = —0.20 Ryd, 
e^f = —0.67 Ryd, = —0.53 Ryd. The lower s valence band should stem from the 2s states of 
N, for which e^ s = —1.35 Ryd; most interestingly, though, the atomic d-states of Ga are found 
almost in resonance with the N 2s states, at ej^f = —1.47 Ryd. From the atomic eigenvalues, one 
can thus infer that the appearance of (i-bands must be expected near the lower valence band. 
Some hybridization of the Ga d states with all other valence states is also to be expected. This 
suggests that the proper treatment of the Ga d states as bands is important for understanding 
the properties of GaN. 



IV. RESULTS 



The structural properties of GaN in the ZB structure from our total-energy minimization are 
a = 8.44 bohr (fi = 150.3 bohr 3 ), B = 1.98 Mbar, and E coh = 10.88 eV/cell (free atom spin- 
polarization included). The reported experimental values of the lattice constant of zincblende 
GaN range |E]|| from 8.50 to 8.58 bohr, so that the experiment-theory discrepancy is around 



1%, which is usual for LDA calculations on semiconductors. We assume the experimental 
cohesive energy per atom pair of the wurtzite phase to be a reasonable estimate of the (as yet 
unknown) experimental cohesive energy per cell of cubic GaN : this amounts to 9.0 eV/cell. 
The overestimation of cohesive energies is a known defect of the LDA 

To follow up the suggested special role of the cation 3d states in GaN, the band structures of 
three compounds are compared in Fig.|2|: those of ZnS, GaN, and GaAs. The role of the metal 
d states is strikingly evident in the GaN band structure (central panel), which turns out to be a 
borderline case between the typical III-V compound GaAs (bottom) in which practically inert 
d bands appear well below the s-like bottom valence band, and the typical II- VI compound 
ZnS (top), in which the d bands lie high in the heteropolar gap. In the case of ZnS, in spite of 
the relatively limited dispersion, a large hybridization with other states takes place (e.g., it has 
been shown [[12| for ZnO that Zn 3d character is as large as 30 to 70 % in valence states usually 
having mostly sp character.) 

It is apparent that a proper treatment of the Ga d shell is essential in determining the band 
structure of GaN. As expected, the d bands are strongly hybridized with the bottom s-like 
valence band stemming from N 2s, which results in a large splitting away from the zone center 
(as sd mixing is symmetry- forbidden at V). For example, at the L point the bottom valence 
band (Li in Fig.||]) and the lower branch of the upper valence bands (L^ in Fig.||D both belong 
to the Li representation of the little group of k, C% v ] the d bands split into two basically inert 
L| doublets ( one of these originates from the non-bonding T^, see Figures in Sec.V), and an Li 
singlet which interacts strongly with the other Li bands (see also the discussion in Sec.V). 

In FigfJ, we show a blow-up of the d bands along the A and A lines. Along A the symmetry is 
C*3v : the Ti2 band goes over to a A 3 state and remains two-fold degenerate, while the Tis splits 
into a A 3 doublet and a A x singlet, which interacts with the bottom and top valence bands. 
Along A, the symmetry is C^ v , so that r 15 splits into a A 5 doublet and a Ai singlet, and T 12 
goes into A 2 and Ai singlets. The compatibility between the latter and the singlet stemming 
from r 15 brings about an anticrossing feature with antibonding and non-bonding character being 
exchanged between the upper and lower level. The splitting of the hybridized band and the 
center of mass of the "inert" d bands is 2.03 eV at X and 1.56 eV at L. It is of some interest 



to note that these features are an amplification of a pattern present, although less pronounced, 
in the d bands of ZnS and ZnSe (we calculated these too, but we do not show them here), in 



accordance with previous all-electron results [13 



The (DFT-LDA) fundamental gap is found to be direct at T, E g = 2.0 eV. The experimental 
value of the fundamental gap is 3.55 eV in wurtzite GaN, and reported values for cubic GaN 
range from 3.25 to 3.5 eV An instant estimate of the self-energy correction to the gap 



plJl^l is provided by a recently proposed simple model [pL5J| , which expresses it as A ~ 9/e eV. 
Since the high-frequency dielectric constant e is unknown as yet for the cubic phase, one can 
use the wurtzite-phase e of 5.8, getting a corrected gap of 3.55 eV. The calculated valence band 
width is 15.99 eV, which is substantially larger than that of the related compounds GaAs and 
GaP (13.0 and 13.1 eV respectively). The main heteropolar gap (s to sp band) is 9.01 eV, while 
the d-to-sp gap is 5.39 eV (the heteropolar gap in GaAs is 3.44 eV). Apart from d hybridization, 
these large gaps are to be expected in view of the ionicity of GaN. 

Several pseudopotential calculations have been recently carried out for GaN, but to our knowl- 
edge there were no all-electron calculation available on it to date for the zincblende structure 
|[(|. A proper comparison of our results with previous theoretical investigations is thus not 
straightforward, due to the problems encountered by pseudopotential methods in treating lo- 
calized nodeless wavefunctions such as those of the N 2p and Ga 3c? states. In the following we 
will compare our calculation with results obtained by others for the zincblende structure only. 



In the mixed-basis calculation of Ref. [0 a pseudopotential including non-linear core correc- 
tions (NLCC) [pl| for the frozen 3c? shell was used for Ga and the Bachelet-Hamann-Schluter 
(BHS) potential |19| was employed for nitrogen. A cell volume of 134 bohr 3 (lattice constant 



a = 8.12 bohr, -3.7 % with respect to our result), a bulk modulus of 2.4 MBar, and a direct 
gap of 2.8 eV were reported. Assuming that the localized part of the mixed basis can describe 
reasonably well the nitrogen valence states despite the low plane-wave cutoff of 14 Ryd, this 
severe underestimate of the cell volume could well be due to the approximate treatment of the 
c?-shell of Ga. This is consistent with past experience with II-VI's [EDI; we return to this question 



in the next Section. We should mention though that the reported cohesive energy is only 8.2 



eV/cell; as has been argued in Ref. |20fl , among the effects of partially active d shells there is 



an increase in lattice constant and a reduction of the cohesive energy (see below), so we would 
expect the pseudopotential cohesive energy to be larger than ours. One may suppose then that 
this calculation may not be fully converged. 

Recently, Palummo et al. [[H|] have studied GaN with plane waves and pseudopotentials, 
treating the Ga d shell as core states, and linearizing the core-valence exchange-correlation 
functional (no NLCC). Their results for zincblende GaN (a tour de force requiring a cutoff of 
120 Ryd) are somewhat puzzling. Using the nitrogen pseudopotential by Gonze, Stumpf, and 
Schemer (GSS) |2| they obtain a = 8.41 bohr (-0.5% from our result), B = 1.69 MBar, and 
-Ecoh = 10.25 eV (-5.7 % from ours), whereas they find ao = 8.15 (-3.4% from our result) and 
B = 2.4 MBar when using the nitrogen BHS potential, in close agreement with Ref. [fT^] . The 



s and p channel of GSS potential and the s, p, and d channels of BHS potentials were used 
respectively, the last channel listed being included in the local potential. As the authors point 
out, this important test signals high sensitivity of the results to the specific potential used for 
nitrogen, due to the high ionicity of the system. The fundamental gap reported is 2.70 eV, 
about 0.7 eV larger than ours [031. This discrepancy could also be attributed to the neglect of 



the Ga d shell (see Sec. IV). 

A further comparison is possible with our calculated frequency of the TO frozen phonon at 
the zone center, which is fi^Q = 600 cm -1 . The experimental frequency is not known in cubic 
GaN, but the analogous mode in the wurtzite structure has ti^o = 533 cm -1 . Although the 
discrepancy is rather large, part of the problem may lie in wurtzite-ZB differences. As an 
example, the TO phonon frequency in cubic SiC is known experimentally to be about 10 % 
higher than in hexagonal SiC [0,0, so this may account for most of the present deviation. On 



the other hand, the pseudopotential result for the wurtzite TO frequency in Ref. |H] suffers 
from a considerably worse discrepancy (fixo = 644 cm -1 ). 

V. EFFECTS OF D ELECTRONS 

The previous Section has demonstrated that the Ga d states play an important role in the 
band structure of GaN, due to their strong hybridization with the N 2s states near the bottom 
of the valence band. Not yet clear is the importance of these active d bands for the cohesive 



properties and the charge density. We consider these questions in this Section, coming to the 
conclusion that GaN is quite similar to the II- VI compounds in these respects. 



A. Cohesive properties 

The question to be answered is whether a reasonable description of the bonding in GaN can 
be obtained with an approach which only treats the sp electrons of Ga and N explicitly as band 
states. In the pseudopotential framework, for example, this could be done by simply treating the 
cation d states as core states and linearizing core-valence exchange-correlation. However, this 
approach has been shown to fail manifestly for the II- VI compounds, giving lattice constants 



which are typically 10 to 15 % too small P0| , ^5| . In this respect, improved pseudopotential 



techniques adopt the so-called non-linear core corrections (NLCC) ||18|| . Hereby the free-atom 
core density is added to the valence density when the exchange-correlation energy is calculated. 
This improves things considerably for the II-VI's, since the lattice constant is now only 3 to 
4.5% below the experimental value [25]. 



More generally, the core states can affect the bonding in three distinct ways. First, if the core is 
frozen to its free-atom shape, an (usually small) error is made which, according to the variational 
principle, must increase the total energy, and hence decrease the cohesive energy. This effect 
disappears when the lattice constant approaches infinity, and it increases monotonically for 
decreasing lattice constant as long as the core states remain very localized on the typical length 
scale of the lattice : in the latter regime this effect leads then to an increase of the lattice 
constant. Second, if the core overlaps appreciably with the on-site valence states, the effect of 
the exchange-correlation non-linearity sets in, resulting in an interatomic repulsion. Third, if 
the core states are large enough to overlap to some extent (or to interact with neighboring-site 
filled valence states, as is the case for the Ga ds and N s in GaN), the closed-shell repulsion 
becomes noticeable. This is a consequence of the increase in the kinetic energy when the cores 
on neighboring sites are made orthogonal. If this contribution is neglected, the lattice constant 
comes out too small and the cohesive energy too large. The magnitude of this effect is governed 
not only by the overlap of the core wavef unctions, but also by the amount of resonance between 
the on-site core eigenvalues for the neighboring sites. This means that the core states can only 



"see" each other if they are in the same energy range. For a system with neighboring atoms of 
different types, it is therefore conceivable that there is a reduced closed-shell repulsion if core 
eigenstates do not resonate. Conversely, for the case of GaN we expect a strong such repulsion 
since, as we have seen, there is a very strong resonance between the Ga 3d and N 2s states. 

To distinguish between the various core contributions to the bonding, we have done additional 
calculations as follows. In each iteration, the core density is obtained by overlapping the frozen 
free-atom cores, and the sum of the free-atom core kinetic energies is taken into the total energy. 
The overlapped core can extend into the interstitial region and into neighboring atomic spheres. 
In cases where the core states are localized enough not to overlap, this "frozen-overlapped- 
core" approximation (FOCA) is equivalent to the usual frozen-core approximation. For GaN, 
however, the Ga 3c? states are so large compared to the interatomic distance that a straight 
frozen-core calculation runs into problems. Of the three core contributions discussed above, 
the FOCA procedure includes only the non-linear exchange-correlation contribution; it does 
not include the core relaxation and closed-shell repulsion effects. It should therefore be closely 
similar to the pseudopotential-plus-NLCC approach, and should shed light on the applicability 
of the latter to GaN and related systems. 

Table I shows the cohesive properties of GaN, ZnS, GaAs, and Si calculated by the all-electron 
FP-LMTO method, both with the full treatment and using the FOCA, and compares them to 
NLCC-pseudopotential results where these are available. For Si, the core is very small so that 
the only neglected term in the FOCA is the core relaxation. In agreement with the discussion 
above, this leads to a slightly larger lattice constant and a smaller cohesive energy. A similar 
behavior is seen for GaAs, for which the Ga d states are energetically well below the As s states 
and, in addition, the lattice constant is rather large. The decoupling of the d states is not 
complete, though, and the relatively small effects of core relaxation and closed-shell repulsion 
are in competition. 

For ZnS as well as GaN, on the other hand, the relevant effect neglected in the FOCA is the 
closed-shell repulsion with the associated underestimate of the lattice constant and overestimate 
of the cohesive energy. Although acceptable values are obtained, the lattice constant is still too 
small by some percent. This is in fact similar to the deviations found when the pseudopoten- 

n 



tial technique with NLCC is applied, as can be seen from Table I. We note though that the 
NLCC-PP errors tend to be somewhat larger. This may be attributed in part to the additional 
approximations needed in the implementation of the NLCC (e.g. the pseudization of the core 
charge inside some inner core radius). 

On the basis of the above considerations we conclude that in the case of GaN no pseudopo- 
tential calculations can be considered reliable if performed without the d shell in the valence 
or (as a minimal requirement) use of the NLCC. This is also the case for II- VI materials, as 
shown above and in earlier work p0| , p5ip6| . Satisfactory results obtained in this context with 



the direct, uncorrected pseudopotential procedure should thus be regarded as incidental pi| .[28 
and possibly due to accidental error cancellation. 

We briefly consider the influence of Ga G?-shell freezing on the band structure of GaN, in 
particular on the direct gap, which is 2.00 eV in the full calculation. Although the absence of 
repulsion between d and top valence p states at zone center causes an increase of the gap (see 



the detailed discussions in Ref. |13|,^(|), a quantitative estimate of gap differences is difficult, 
due to volume dependences; the high bulk modulus and deformation potentials of the relevant 
states cause effects similar to those observed for other semiconductors ]2"9fl. A gap of 2.8 eV is 



reported in Ref. |T7j|, and one of 2.7 eV in Ref. [2~I|. In our frozen-core calculation we find a gap 



of 2.66 eV at the frozen-core theoretical lattice constant (8.30 bohr). When calculated at the 
theoretical lattice constant of the full calculation (8.44 bohr), the frozen-core gap is 2.20 eV, 
10 % larger than the full gap at the same volume. We report in passing our estimate for the 
deformation potential of the lowest conduction state, 5 meV/kBar, in good agreement with the 



experimental value of 4.7 meV/kBar |]T6 



As a summary of this Subsection, the core states can influence the bonding properties in 
three ways, namely through non-linear core- valence exchange-correlation effects, core relaxation, 
and closed-shell repulsion. To obtain a resonable description of bonding in II- VI and III-V 
semiconductors with d states in the valence range, it is imperative to include non-linear core 
corrections in a pseudopotential treatment. The core relaxation has a minor influence and can 
normally be neglected. The closed-shell repulsion is rather significant in those cases for which 
the cation d states are in the same energy range as the anion sp states. In this sense, GaN turns 



out to be very similar to the II- VI compounds. For all of these materials, a full treatment of 
the d electrons as band states is needed in order to reduce the remaining 3-4 % underestimate 
of the lattice constant to the LDA-standard 0.5-1%. 



B. Wavefunctions and densities 

The similarity of GaN and the prototypical II- VI ZnS has been observed in the band structure 
and cohesive properties. Such effects are due to the cation d shell being not inert in these 
systems. To render this visually, we present a selection of partial charge densities for GaN. 
This is similar in spirit to the work of Wei and Zunger |p0[ , who have investigated in some 
detail the effect of cation d states on (among others) the charge densities of selected states 
and on the total densities in several II- VI compounds. We will be considering charge densities 
at specific k-points (the modulus squared of the wave function ipi (k, r) at k) as well as k- 
integrated band contributions to the density. All charge densities are plotted in a (110) plane. 
For convenience, superscripts d and v are used to label states originating from d and top valence 
bands respectively, whereas states with no superscripts are from the bottom valence band (see 
Figfjfor a guide to the relevant states). Density plots have lowest contour and contour spacing 
of 0.5 and 1 electron/cell, respectively. 

In the zincblende structure, the mixing of p and d states is permitted |20| even at the T point, 



which has Td as group of k (contrary to the situation in the diamond structure, where it has 
the full Oh symmetry). The upper valence bands at T will no longer have pure p character, d 
character being admixed into them and, conversely, the (i-bands will contain p components. The 
mixing of s and d states is not allowed at T, but at a general point of lower symmetry, states with 
all symmetries of the rotation group will in general be allowed to interact. Therefore, charge 
density and density of states (which are obtained by averaging over the Brillouin zone), and the 
full band structure will be only partially reminiscent of the specific symmetry of the original 
atomic states. This is the case, e.g., for the Tf 5 triplet of d origin which is only p— admixed 
at the zone center, but gives rise to a strongly s— and p— admixed singlet away from the zone 
center, and for the Ti s-like singlet, which is similarly admixed with pd character away from 
T. Already at the zone center, the ordering of the — like r 15 triplet and T 12 doublet (which 



remains practically inert throughout the zone, resulting in a non-bonding state) is reversed with 
respect to what is expected 0] for a tetrahedral environment, due to pd repulsion with the top 
valence band (qualitative schemes of this interaction are given in Ref. ||20| ). 

In Fig.[|we plot the squared wavefunctions at the zone center for the states rf 5 ,rf 2 and r^ 5 
in GaN. The former two states originate from the tetrahedral splitting of atomic d states (the 
5-fold degenerate d multiplet is split in a Ti 5 triplet and a r i2 doublet), the latter is the bonding 
p— like valence state. The Ti2 state remains very similar to a c^-^-like atomic orbital and forms 
a non-bonding, inert state. In the other states, the prevailingly ionic character of the bonding 
is visible; the admixture of d xz -like and p character results in a weak bonding contribution to 
r^ 5 and an antibonding one to r" 5 . 

While at zone center only pd mixing is allowed, this restriction does not apply to lower- 
symmetry points. We consider as representative the squared wavefunctions at the L point. At 
that k-point the bottom valence band is of Li symmetry; the Ij 5 triplet splits further into an L3 
doublet, analogous to the (r 12 )-l4 and an Lf singlet. Figure |5] shows the symmetrized squared 
wavefunction at L (the average over the eight equivalent L-points in the Brillouin zone), for the 
states Li, L^, L", L3 (top to bottom). We notice that the bottom valence Li has a bonding d 
admixing, and Lf and L3 are instead antibonding. 

Next, in Fig.|6| we plot the k-integrated partial charge densities nj(r) for the bottom, d, and 
top valence bands. The first two are sd-\ike and have some bonding contribution from the 
d-like states, the latter is ionic p-like and has an antibonding <i-admixing, which is visible as 
a depletion along the bond direction. As expected the Brillouin zone summation has brought 
about a mixing of all angular momenta (in particular sd in the two lower bands). The top 
valence states are clearly dominated by ionic charge transfer. 

Finally, we compare GaN to ZnS and GaAs in Fig.[7[ Shown for the three materials (top 
to bottom : ZnS, GaN, GaAs) are the states r^ 5 (right) and the k-integrated densities of the 
bottom valence band (center) and of the top valence band complex (left). The similarity of GaN 
and ZnS is remarkable. The effect of d admixing is more visible in the ZnS top valence charge, 



as the ci-bands are close to and significantly hybridzed with the top valence bands The 



bonding contribution of the d electrons to the ionic bottom valence charge is instead stronger 



in GaN, as here the hybridization of s and d bands is larger. The top valence band, zone- 
center state ry 5 (right) shows a conspicuous antibonding contribution of the d electrons, of a 
similar kind in ZnS and GaN. Even in GaAs, though, a slight <i-admixing is visible in this state 



(analogous to that observed recently in InP, InAs, and InSb ||31||), and its bottom valence charge 
also shows some d component. We should mention that ionicity is only partially responsible 
for the differences between GaAs and GaN, or the similarities of ZnS and GaN, as the ionicity 
of GaN is 0.43 as compared to 0.33 of GaAs and 0.70 of ZnS [fl6 1 . The features of the partial 



densities should also be directly compared with the results of Ref. |2(| for II- VI compounds. 

We conclude that in contrast to most other III-V compounds (represented here by GaAs), 
GaN does exhibit features of the charge densities and band structure which resemble those 
observed [^0J for II-VI's (represented here by ZnS). This is due to the d shell of Ga being not 



inert in this compound, as a consequence of the energy resonance of N s and Ga d levels. 

C. Density of states 

Although mixing of all angular momenta takes place due to symmetry lowering in the crystal 
environment, the content of a specific angular momentum character of a crystal state may 
be recognized in several ways. As seen above, inspection of partial charge densities is useful 
in this respect. Another technique allowing a monitoring of the site and angular momentum 



dependence of the charge distribution is the Mullikan decomposition |32| of the total density of 
states, which is natural to the LCAO-like basis used here. The total DOS can be decomposed 
into site and angular momentum contributions, and if desired, into bonding and antibonding 
overlaps. 

The LMTO basis set naturally gives a decomposition of the calculated wavefunctions as a sum 
over sites and angular momenta as ipk = J2 U l V'iM j where v is a site index and L a composite 
angular momentum index. We then take the decomposition of the norm as weight to add into 
the partial densities of states as 

i = ^ki^k) = E Re [^kiri1)]. (i) 

v,L 

The advantage of this approach over the common technique of projecting onto angular momenta 



inside the spheres is obvious : first, the decomposition is essentially independent of all sphere 
radii, and second, no DOS contribution is associated with the empty spheres or the interstitial 
region. 

The atom-projected DOS for s, p, and d angular momenta are given in Fig.[| A gaussian 
broadening of 0.02 Ryd and a mesh of 344 special points in the irreducible wedge of the Brillouin 
zone are used. The s and d bands have contributions from both Ga and N sites, but this is a 
consequence of their having exchanged characters to a considerable extent. Most s character 
comes from N, whereas most of the d character resides on the Ga site, and is barely sizable on 
N over the whole valence band span (the total charge of d character contributed by the N site 
is about 0.08 electrons). The N p states dominate the top valence band, although of course spd 
Ga contributions are visible. The s DOS shows a predominant contribution of N 2s to both of 
the lower bands, and almost none in the upper bands, where some Ga s character is present 
instead. 

VI. SUMMARY 

We have calculated the properties of /5-GaN using the ab initio all-electron full-potential 
LMTO method. Agreement is found with existing experimental data, and predictions are made 
for other quantities. The explicit consideration of the d electrons of Ga has been shown to be 
of importance in determining the properties of the material. Due the resonance of N s and Ga 
d states, the band structure, DOS, and charge density of GaN are quite peculiar for a III-V 
compound, as significant 3d character is present in all the valence bands. Several features of 
GaN of the band structure, charge density, and density of states do actually resemble those of 
II- VI materials. 

A frozen-overlapped-core procedure was used to identify the way in which the core and d 
states influence the cohesive properties of GaN as compared to some other relevant materials. 
We deduce that, although the non-linear exchange-correlation correction is the most important 
effect in GaN and ZnS, the closed-shell repulsion also plays a significant role; neglect of this 
effect leads to a lattice constant which is too small by some 3% in both systems. This sets an 
upper limit on the accuracy which can be expected of methods which do not treat the d states 



esplicitly. Therefore, even pseudopotentials including non-linear core corrections, although going 
a long way in the right direction, are generally not sufficient for a description of GaN, as well 
as of ZnS and similar system, at the standard LDA level of accuracy for semiconductors. 
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FIGURES 

FIG. 1. Non relativistic free- atom LDA eigenvalues for some elements involved in the formation of 
III-V and II- VI compounds. 

FIG. 2. FP-LMTO band structure of GaN (center), compared to those of ZnS (top), and GaAs 
(bottom). State labels for the central panel apply straightforwardly to the others. The energy zero is 
the valence band top. The bands are calculated at the theoretical lattice constants. 

FIG. 3. Detail of the d bands of GaN along L-F-X. Degeneracies are indicated in parentheses. 

FIG. 4. Squared wavefunctions of the states rf 5 (top), Tf 2 (center), rf 5 (bottom) at the zone center 
in GaN. 

FIG. 5. Symmetrized squared wavefunctions of the states Lg, L\, Lf, and Li (top to bottom) at 
the L point in GaN. 

FIG. 6. Partial k-integrated charge densities in GaN for the top valence bands (top), d bands 
(center) and bottom valence band (bottom). 

FIG. 7. Comparison of ZnS (top row), GaN (center) and GaAs (bottom). Right column : state 
F\ 5 ; central column : bottom valence band partial density; left column : top valence band partial 
density. 

FIG. 8. Mullikan angular momentum decomposition of atom-projected densities of states. From 
top to bottom : Ga s, p and d states, and N s, p and d states. 
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TABLES 







all-electron a 


FOCA 


b 


NLCC-PP C 


Exp. d 


GaN 


a 


8.44 


(-1.2) 


8.30 


(-2.8) 


8.12 


(-4.7) 


8.54 




B 


1.98 




2.00 




2.40 




- 










0.33 










ZnS 


a 


10.13 


(-0.8) 


9.92 


(-2.9) 


9.80 


(-4.0) 


10.21 




B 


0.75 




0.83 




1.06 




0.77 




<$-Ecoh 






0.34 










GaAs 


a 


10.62 


(-0.6) 


10.71 


(+0.2) 


10.50 


(-1.7) 


10.68 




B 


0.75 




0.65 




0.77 




0.77 










0.08 










Si 


a 


10.22 


(-0.4) 


10.24 


(-0.2) 






10.26 




B 


0.96 




0.96 








0.99 




5E coh 






-0.41 











TABLE I. Calculated structural properties for GaN, ZnS, GaAs, and Si. The lattice constant 
ao is in atomic units, B is the bulk modulus in Mbar, 5E co h = E^^ en — -E^ol! ^ s the difference of 
cohesive energies (eV/cell) in the frozen-core and full all-electron calculations, respectively. Numbers 
in parentheses indicate deviations from experimental data. 
a : this work, full FP-LMTO 
b : this work, frozen-overlapped-core FP-LMTO 



c : pseudopotentials with non- linear core corrections (Ref. |l7| for GaN, Ref. 25], for ZnS, Ref. 
for GaAs) 

d : data from Landblt-Bornstein : Numerical Data and Functional Relationships in Science and Tech- 
nology, New Series Group III, K.- H. Hellwege and O. Madelung eds., (Springer, New York, 1982), 
Vol. 17a and 22a and references therein. For cubic GaN an average of the reported experimental data 
J2| is used. 



